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Abstract: In this paper we present the results of the analytical and numerical studies of 
the plume interaction with the base flow in the presence of base out-gassing. The 
physics-based analysis and CFD modeling of the base heating for single solid rocket 
motor performed in this research addressed the following questions: what are the key 
factors making base flow so different from that in the Shuttle [1]; why CFD analysis of 
this problem reveals small plume recirculation; what major factors influence base 
temperature; and why overheating was initiated at a given time in the flight. To answer 
these questions topological analysis of the base flow was performed and Korst theory 
was used to estimate relative contributions of radiation, plume recirculation, and 
chemically reactive out-gassing to the base heating. It was shown that base bleeding 
and small base volume are the key factors contributing to the overheating, while plume 
recirculation is effectively suppressed by asymmetric configuration of the flow formed 
earlier in the flight. These findings are further verified using CFD simulations that 
include multi-species gas environment both in the plume and in the base. Solid 
particles in the exhaust plume (A1 2 0 3 ) and char particles in the base bleeding were also 
included into the simulations and their relative contributions into the base temperature 
rise were estimated. The results of simulations are in good agreement with the 
temperature and pressure in the base measured during the test. 
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1 Introduction 

Recent test flight with single motor first stage design posed significant challenge in predicting base 
environment using both shuttle space flight database and CFD code. The post flight detached LES 
performed for this configuration revealed plume recirculation at the level substantially below one required 
for observed overheating. To gain further insight into root causes of the observed overheating and to 
facilitate the prediction of the thermal environment for the next generation of the heavy lift vehicle the 
following analysis was performed. 

2 Problem Statement 

A model of the vertically stacked rocket with single motor first stage was built in CFX and simulated for 
a number of time instants during the flight. The base flow sub-domain was next simulated with increased 
resolution. Mesh adaptation was used to resolve flow mixing in free share layers as shown in Fig. la. The 
k-co based Shear-Stress-Transport (SST) model and BSL Reynolds stress model were used in simulations. 
The results of the simulations in the absence of the base bleeding and A1 2 0 3 particles in the plume were in 
agreement with earlier detached LES performed for this geometry by one of the authors (FC) showing 
base temperatures corresponding to the aerodynamic heating and plume recirculation at the level of few 


percent. Further insight into the root causes of the observed phenomenon was provided using topological 
analysis of the base flow that demonstrated strong asymmetry of the flow, when the region of plume 
recirculation with the relatively high temperature is very small and pressed against the plume boundary 
preventing plume recirculation mixing with the main volume of the base. This configuration of the 
trapped based flow is formed earlier in the flight. Next, relative contributions of the radiation, plume 
recirculation, and base bleeding into the base heating were estimated using Korst theory [2], showing that 
base bleeding may be one of the key factors that determine base overheating observed in the test flight. 



Figure 1: (a) Mesh Adaptation showing the location of mixing layers; (b) Char and A1 2 0 3 particle 
trajectories are shown together with the contour plot of the temperature in the plane of symmetry. 


To verify these findings a model of chemical kinetics of the pyrolysis of the TPS was developed [3] to 
determine chemical composition, temperature, and mass flow rate of the base bleeding. Parameters of the 
bleeding determined from the chemical kinetics analysis were next substituted into the CFD simulations. 
The results of the simulations of the base flow with base bleeding confirmed significant increase of the 
base temperature in agreement with the test flight data. 


2 Conclusion and Future Work 


Theoretical analysis and CFD simulations performed in this research revealed that the base heating 
anomalies observed in a recent test flight can be attributed to the following factors: small volume of the 
base flow in a single motor first stage configuration; significant base bleeding of hot chemically active 
components and char trapped in the base; and strong asymmetry of the base flow with relatively small 
region of plume recirculation. An additional important factor that influence the overheating measured 
during the flight is a strongly non-uniform distribution of the hot flow in the base, which is mainly 
concentrated along the base walls for both base bleeding and plume recirculation. 

The methodology of the base flow analysis developed in this research is currently being transferred, 
extended, and adopted for base flow in nozzle clusters. 
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1 Introduction 

The launch vehicle and missiles encounter base heating due to the interaction of free stream and the 
rocket exhaust. Current research attacks this problem on several directions. Modern experimental 
techniques [l]-[8] allow for visualization of the base flow with nanosecond resolution [8] including 
visualization of large coherent turbulent structures [3], [4]. The analysis of small scale experiments and 
post-flight databases underlies recent progress in the development of the numerical methods [12], [13] 
suitable for more accurate base flow predictions. The progress in the CFD analysis of the base flow is 
paralleled with the development of the numerical methods for analysis of the plume radiation. At the 
same time recent years witnessed renewed interest to the theoretical analysis [16]-[19] of the base flow 
interaction with rocket plume that continues the line of reasoning developed in the Chapman-Korst theory 
[20]-[23]. 

Despite this progress, accurate predictions of the base temperature and heat fluxes during the ascent 
remain a challenge to computational fluid dynamics. Difficulties stem from the fact that the fundamental 


flow multi-physics at the base is not well understood. Often the interplay of the base flow with the plume 
recirculation, afterburning, chemically active outgassing, and radiation has to be taken into account 
rendering the problem of the base flow calculations a formidable task. 

In particular, recent test flight with single motor first stage design posed a significant challenge in 
predicting base environment using both semi-empirical approach relying on the shuttle space flight 
database [23] and CFD solution. Earlier extensive numerical analysis including detached LES simulations 
revealed unexpectedly that the level of the plume recirculation is substantially below the one required for 
explanation of the observed overheating. At the same time post flight data analysis indicates significant 
outgassing at the base that could substantially contribute to the observed increase of the base temperature. 
In this paper we present results of the research performed to analyze the root causes of the observed 
overheating and to facilitate the prediction of the thermal environment for the next generation of the 
heavy lift vehicle. 

The paper is organized as follows. In the next section the problem formulation is presented, the mesh 
generation, including brief description of the computational methodology, mesh generation and 
adaptation. The results of the analysis of the base flow with neglected radiation and the base outgassing 
are presented in Sec. 3. The comparison of the base flow topology observed in numerical simulations with 
the predictions obtained using recent advances in the Korst theory is given in Sec. 4. The model of the 
base outgassing and the base bleeding characteristics are discusses in See. 5. The effect of outgassing and 
afterburning on base environment are discussed in Sec. 6. Finally, a brief summary of the obtained results 
and discussion of the future work is given in Conclusions. 


2 Problem Statement 

In a recent test flight of a vertically stacked rocket with single motor first stage a significant overheating 
of the base environment has been observed. The base temperatures predicted using both CFD analysis 
[10], [11] and semi-empirical approach [23] significantly underestimate observed base temperature. A 
deviation of the predicted temperatures from the measured data is found already earlier in the flight. At 
some points of the flight-trajectory absolute deviations are as high as 65-85% of the predicted values. 
Similarly, CFD code significantly underestimates observed convective heat transfer coefficient, while 
semi-empirical methodology provides accurate predictions of the heat transfer coefficient for the second 
half of the flight. At the same time radiative heat flux predicted using reversed Monte-Carlo code agrees 
well with the observations in the first half of the flight. 

An extensive hybrid RANS/LES simulations performed with substantially improved mesh resolution of 
the base flow confirmed earlier predictions for the base temperature obtained using Member’s baseline 
model. One of the most important results of the RANS/LES analysis was an observation of relatively 
small plume recirculation in the base. The predicted levels of the plume recirculation were up to 5 times 
smaller than levels required for explanation of the observed overheating. 

These observations posed significant challenge to the development of the predicting methodologies. 

The research presented in this paper addresses the following questions: what are the key factors making 
base flow in the single motor first stage so different from that in the Shuttle; why CFD analysis of this 
problem reveals small plume recirculation; what major factors influence base temperature; and why 
overheating was initiated earlier in the flight. 

An important clue to the possible explanation of the measured overheating was an observation of the 
strong outgassing of the thermal protection system (TPS) is clearly seen on the flight video. Therefore a 
significant part of the research was devoted to the analysis of the thermal decomposition of the TPS and 
the changes it introduces in the base flow. We now describe in more details computational methodology 
and mesh properties of the model. 



3 Computational Methodology 


Commercial CFX and FLUENT software capable of solving diverse and complex multi-dimensional fluid 
flow problems was used in this research. Governing equations for the turbulent compressible gas flow 
were taken in the form 
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Note that the last term in this equation that involves the divergence of velocity 2/ 3 /LL eff dU k /dx k is 


neglected in ANSYS CFX. 

One of the key features of the base flow that has to be captured in simulations is the flow separation. 
Accordingly SST model of Menter was chosen as a baseline model following recommendations of the 
NASA Technical Memorandum [30]. Advanced formulation of the SST model is currently available in 
the turbulence modeling program lead by Dr. Menter at ANSYS CFX [31]. In this formulation the 
standard Wilcox model is blended with transformed k-s model resulting in a well-known set of equations 
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with the following set of parameters /? =0.09, ai- 5/9, /?i=0.075, o k i- 2, o m i-2, a 2 =0A4 , /? 2 = 0.0828, a k2 = 1, 
^2=1/0.856. The transport behavior is obtained by a limiter to the formulation of the eddy-viscosity = 
paiMnax(a7CU,5F 2 ) and a standard choice of the blending functions Fj and F 2 . 

The CFX wall boundary treatment exploits the robust near-wall formulation of the k-co model and 
switches automatically from a low-Reynolds number formulation to a scalable wall functions treatment 
based on grid density. This so-called scalable wall function treatment of the SST model gives more 
accurate results for a wider range of grid densities. 

The CFX solver combines the geometric flexibility of finite element methods with the conservation 
properties of the finite volume by integrating conservation equations over a control volume. The pressure 
integral terms in the momentum integral equation and the spatial derivative terms in the integral equations 
are evaluated using the finite-element approach. An element is described with eight neighboring nodes. 
The method is using upwind differencing to ensure global convergence. 

An important feature of the ANSYS solver is the ability of adaptive refinement of the unstructured mesh 
using more than one solution variable. For example, the adaptation criteria, A u for a given mesh edge i of 
length /,, is of the form 


4 = 1 , 




where fa is the j-th adaptation variable Afa is the global range of this variable over all notes, Afa t is the 
variation of the adaptation variable on a given edge, and N ^ is a scalar for adaptation variable. For more 



than one solution variable, the adaptation criteria are calculated for each variable at each edge, and the 
sum over all adaptation variables is used. 

It is essential in the context of the base flow research that ANSYS solvers allow modeling of multi- 
component flow with arbitrary number of species. A number of options is available in ANSYS for 
modeling chemical reactions. However, the rate of chemical reactions in the base is limited mainly by the 
turbulent mixing. Accordingly, eddy dissipation model and eddy dissipation combined with finite-rate 
chemistry were used to model afterburning of the TPS outgassing in the base. 

The solid phase was modeled using Lagrangian tracking of a number of particles through the fluid for 
both Al particles in the plume and char particles in the base flow. For the drag coefficient a limiter is 
introduced [28] to ensure correct behavior in the inertial regime that has the form 
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The particles temperature is determined by the convective heat transfer, chemical reactions, and the 
absorption/emission of radiation at the particle surface 
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The oxidation of the char particles in this model is coupled with eddy dissipation calculation for the 
combustion of the volatile gases in the gas phase. 

Overall ANSYS solvers offer a unique combination of capabilities suitable for the analysis of multi- 
physics phenomena determining the base environment using limited resources. 

The analysis of relative contribution of various heat sources to the base temperature was performed in 
several steps. Initially the computation was performed for the whole computational domain in a number 
of points along flight trajectory. The chemical species transport with frozen chemistry equations were 
solved for both air free stream and plume flow. The chemical composition and distribution of the species 
(including distribution in the boundary layer near nozzle wall) in the plume flow were provided by the 
standard NASA code for plume flow field and specified downstream of the nozzle throat at a plane cut 
through the nozzle flow approximately 1 m away from the nozzle exit. 

At the next step the computational domain was reduced to the analysis of the plume free stream 
interaction at the nozzle exit. The Al particles transport coupled by the momentum and heat transfer to the 
gas plume flow was added to the model. 

Analysis of the convergence of the predictions of the pressure, temperature, and the species distribution as 
a function of mesh size was performed at the next step. The main attention at this step was paid to the 
resolution the free shear layers separating from the racket aft skirt and from the nozzle edge and the 
recompression region including recompression shock. It was important at this step to implement the 
incremental adaptation capabilities of ANSYS CFX known as hierarchical refinement of unstructured 
mesh. 

Once convergence of the base temperature and base pressure predictions was achieved an extensive 
analysis of the base for structure, base bleeding, and afterburning was performed to explain low levels of 
plume recirculation and high base temperatures observed in the test flight. 


4 Model Development 

4.1 Geometry and Grid Generation 

Computational domains and a typical grid structure for a model of the vertically stacked rocket with 
single motor first stage are shown in Figure 1 . Initially unstructured grid with size of a few million cells 



was generated using CFX for the whole computational domain. The flow field around the rocket was 
calculated for several points along the flight trajectory. Next, to improve the resolution of the base flow 
the computational domain was restricted to the region near the nozzle exit where free stream interacts 
with plume as show in the Figure 1 (right). The inlet boundary conditions were provided by the solution 
of the problem in the whole computational domain at the previous step. The fully coupled (except 
radiation) transport equations for the Al particles were added to the model at this stage. 



Figure 1 : Computational domain for the whole rocket is shown on the left. The computational domain 
restricted to the flow field around the nozzle is shown on the right. 


4.2 Near wall flow resolution 

To resolve boundary flow in the near wall regions automatic near- wall treatment for co - based models was 
used [28]. According to this approach the boundary layer should be resolved with at least 10 mesh points 
within the layer for scalable wall functions. The following correlation is used in ANSYS to estimate the 
width of the boundary layer { S) in a given computational domain with characteristic length L and 
Reynolds number Re 

8 oc0.035LRe" 1/7 . 

The estimated width of the boundary layer for both nozzle flow and base flow was of the order of 0.01m. 
To satisfy the criteria of the minimum number of nodes for wall function the initial thickness of the mesh 
layer was taken 5xl0" 4 m and at least 10 layer where placed within boundary layer, while the total number 
of the boundary layers was 30, as shown in the Figure 2. No attempt was made to resolve heat transfer 
coefficient at the nozzle wall at this stage of the analysis. 



Figure 2: (left) Boundary layers at the nozzle surface, (right) The y+ values at the nozzle surface. 




The low-Reynolds number implementation of the model could not be achieved for all the values of used 
parameters. And automatic near wall treatment that allows for a smooth shift from a low-Reynolds 
number form to a scalable wall-function formulation was used as default throughout this research. 


4.3 Mesh adaptation 

In order to resolve free shear layer near the nozzle exit hierarchical refinement of unstructured mesh Mesh 
adaptation was accomplished in ANSYS CFX. Mesh was selectivity refined in 2 or 3 steps in areas 
marked for adaptation by the solver depending on specified adaptation criteria. The Mach number was 
selected as the adaptation variable to guarantee stable and quickly converging solutions. Mesh density 
was automatically increased in locations where solution variables change rapidly (see Sec. 3 and [28]). 
Mesh adaptation did not add extra levels to the boundary layers, but rather it added and subtracted nodes 
within the plane of each layer. Accordingly, mesh adaptation did not change the model’s y+ resolution. 



Figure 3: Mesh refinement in the base flow region obtained in two steps of the hierarchical refinement using Mach 
number variable in adaptation criteria, (left) Mesh before refinement, (right) Mesh after refinement. 


The results of the two-step mesh refinement area illustrated in Figure 3. The final mesh resolves free shear 
layers both on the side of the free stream and on the plume side. The size of the mesh depends on the 
point along the flight trajectory, and in this particular case it was a few million cells. 


4.4 Transition to axisymmetric model 

The size of the mesh after adaptation for 90° nozzle domain shown in Figure 1 (right) could be as high as 
hundred million of cells. To meet the limited computer resources restrictions and in the view of the fact 
that axial asymmetry of the problem at hand is relatively weak, the computational domain was further 
reduced to the analysis of a 5° sector. We note that the physics of the base flow in this case is primarily 
affected by the resolution of the free shear layers. To resolve these layers a very high local mesh density 
is required near the nozzle edge as shown in the Figure 2 (left). Note also that using axisymmetric quasi- 
3D solution in 5° nozzle sector with selective adaptation based on the Mach number variation allows one 
to reveal important features of the base flow that would otherwise require mesh size of several hundred 
millions cells to archive similar resolution using nonselective mesh refinement in the 90° nozzle sector. 
The results of the analysis of the base flow obtained with this mesh are presented below. In the next 
section the results of the calculations of the base flow without base bleeding is discussed, the base flow 
topology is revealed, and the unexpectedly how levels of the plume recirculation observed in the earlier 
simulations for the vertically stacked rocket with single solid motor are confirmed and explained. The 
discussion of the base bleeding characteristics and the base bleeding effect on the base temperature will 
be presented in Sec. 7. 




5 Results of the simulations without base bleeding 


5.1 Structure of the base flow 


Results of the simulations of the free stream interaction with the rocket jet obtained using the model 
described in the previous section for one point along the flight trajectory are shown in Figure 4 and Figure 
5. A high-contrast contour plot of the Mach number distributions reveals clearly the well-known structure 



Figure 4: Contour plot of the Mach number: (a and a') mixing layers; (b) barrel shock of under- 
expanded jet; (c) recompression shock waves near wake; (d) base flow; (e) plume boundary; 

(f and f ) expansion fans. 

of the flow field corresponding to the free stream interaction with the rocket jet. The two shear layers (a 
and a') bounding the base flow (d) on the plume side and the free stream side can be seen in the figure. On 
passing the rocket skirt the free stream is changing its direction through expansion fan (f) to match the 
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Figure 5 (left) Pressure predictions for various flight trajectories are shown by the open circles in 
comparison with measured data (blue and green solid lines) and total free stream pressure (red line), 
(right) Temperature predictions are shown by the open circles in comparison with measured data (green 
solid lines) and total free stream temperature (blue line). 

base pressure. The free stream boundary layer separates from the skirt and forms free shear layer (a) that 
follows free stream direction. Similarly, on the plume side the under-expanded jet is expanding through 


an expansion fan (f) to match base pressure. The jet boundary layer separates from the nozzle at the 
nozzle exit and forms free shear layer (a') that follows the direction of the plume boundary. The two shear 
layers merge in the recompression zone bounded on both sides by the recompression shocks (c) coupled 
to the shear layers. As a result, the base flow is trapped between two shear layers. 

The results of the calculation of the base pressure for mesh adapted at various points along the flight 
trajectory are shown in the Figure 5 (left). It can be seen from the figure that the predictions for the base 
pressure follow very closely the total free stream pressure in good agreement with earlier calculations 



Figure 6: (left) Temperature distribution, (right) HC1 mass fraction distribution. 


[10] , [11], confirmed by the RANS/LES calculations with high mesh density for one specific trajectory 
point. The pressure predictions are also in good agreement with the measured data for the first part of the 
trajectory, but underestimate measured pressure later in the flight. Furthermore, calculated pressure fails 
to predict the crossover of the base pressure with ambient pressure observed during the flight test. Note, 
that such a crossover usually indicates the mass addition to the base related to the initiation of the plume 
recirculation [25]. 

The predicted temperatures are shown in the Figure 5 (right). It can be seen from the figure (see also 
contour plot of the temperature distribution in the base flow shown in Figure 6 (left)) that predicted 
temperatures follow closely the total free stream temperature in agreement with earlier calculations [10], 

[11] . Both earlier calculations and current predictions significantly underestimate observed base 
temperature already earlier in the flight. At some points of the flight-trajectory absolute deviations are as 
high as 65-85% of the predicted values. 

Earlier research [25]-[27] has shown that the main source of the elevated base temperature is related to 
the plume recirculation. For example, in the Shuttle [25] or in Titan III [27] plume recirculation is 
initiated when plume of the two solid rocket motors overlap. In the single motor first stage design plume 
recirculation may occur only via mass transfer in the mixing layer on the plume side. To explain the 
observed base overheating plume recirculation at the level approximately 10% of the mass fraction of the 
base flow is required. However, our calculations show (see Figure 6 (right)) a much smaller level of 
plume recirculation, which is approximately 1% of the mass fractions the base flow. This result is in 
agreement with the earlier calculations [10], [11], verified by the RANS/LES calculations for one specific 
trajectory point. 

To understand unexpectedly low levels of plume recirculation in this geometry let us consider base flow 
topology in more details. 

5.2 Topology of the base flow 

Topological insight into the key characteristics of the base flow can be obtained using recent advances 
[16]-[19] of the Korst theory [20] -[24]. According to this theory the flowfield within the shear layers is 


divided by two main streamlines r(r ) and s(s ) for the free stream (plume). The r and r ' joint in the end of 
the wake at stagnation point Pj in the recompression region as shown in Figure 7 (left). The streamline s 
terminates in the stagnation point P 2 in the base, while streamline s ' propagates downstream along the 
plume boundary. The importance of the asymmetry of the base flow for the free stream and plume having 
different total pressure was emphasized in [18], [19], where it was noted that, in general, when two free 
jet surfaces join, the stronger flow (higher total pressure) is the ejector, and the weaker is the feeder. 



Figure 7: (left) Sketch of the main streamlines separating various regions of the base flow, (right) vector 
plot of velocity distribution of the base flow colored by the velocity magnitude. Solid lines indicated 
main streamlines. Solid lines indicated main streamlines separating various regions of the base flow. 

This asymmetry is responsible, in particular, for the low levels of the plume recirculation in the current 
geometry of the base with a single rocket motor. The latter can be understood by identifying location of 
the key separating streamlines and stagnation points in the vector plot of the base flow velocity as shown 
in Figure 7 (right). It can be seen from the figure that the pattern of the streamlines r ' and s'i s degenerate 
and strongly pressed against the plume boundary because of the very large difference between the total 
pressures of the free stream and rocket jet. Note that the streamline s impinges the base wall at stagnation 
point P 2 located at the nozzle exit. We see that the outflow in the corridor s '-r 'separates mixing layer on 
the plume side from the base flow preventing any significant plume recirculation in the base. Therefore 
the base flow in a given geometry corresponds to the trapped vortex bounded by the streamline s. 

Further analysis reveals that transition to the trapped flow occurs earlier in the flight, which correlates 
well with the initiation of the base overheating observed in the test flight. 

To estimate contribution of other physical mechanisms to the base heating let us briefly recall the mixing 
properties of the shear layers boundary the base flow. The equations for the turbulent shear layer have 


the form 
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By solving the momentum equation in the incompressible approximation for free mixing layer using 

Korst theory [18] -[24] one obtains 
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where the compressibility corrections are introduced through the experimentally verified correlations 
[24], [24] 
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Temperature and pressure distributions can be obtained from Crocco's integral of the boundary-layer 
equations 
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Once temperature and velocity profiles are established the mass and heat in- and outflows can be 
determined by integration across streamlines in the mixing layer. For example 
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Notice, however, that in the first approximation the flow in the corridor between streamlines s-r , s-r\ and 
s -redoes not enter the region of the dead air corresponding to the trapped vortex in the base. 

It is clear that in the first approximation of the Korst theory there is no ventilation of the air trapped in the 
base of a single motor rocket. Accordingly, the base temperature is prompt to significant deviations from 
the total free stream temperature due to various secondary sources of heat. For the air trapped in the base 
with mass M p and heat capacity c p the energy conservation equation reads 


CpM h ~dt~® bb+ ® r + ® ch - ® Ioss + K ^dx 


(7) 


where the following contributions are taken into account: Q bb is the heat flux due to the base outgassing 
and chemical reactions, Q r is the heat flux due to the radiation, Q ch is the heat flux due to the char flux and 
afterburning. The losses term Q ioss is introduced to take into account mass deposition from the base into 

the mixing layer. The existence of much deposition is clear from the fact that. The term fcdT /3x|^ is the 

heat flux due to heat conduction at the boundary of the trapped vortex. The last term is automatically 
included into the base flow calculations described above. 

Note that substantial amount of the outgassing was visually detected in the base during the flight test. 
Therefore in what follows the model is extended to include contributions of the outgassing Q ch and char 
burning Q bb into the base temperature calculations. 


6 Ablation of thermo-protection system 

To estimate the mass and heat fluxes in the base due to the base bleeding a model of decomposition of the 
thermal protection system was developed. The TPS typically consists of Polyurethane (PU) foam, 
catalysts, blowing agents, surfactants, and flame retardants. In extreme heating environments caused by 
the plume radiation and heat convection the TPS ablates producing hot gas and smog [32]-[35]. The PU 
decomposition is a complex process that includes three main steps. First, endothermal pyrolysis and 
internal oxidation of the PU occurs. The energy absorbed during this process is H . The products of this 

reaction are thermo-char with density p , gas and oil droplets (yellow smoke) with total density p . At 

1 c gp 

the second stage, pyrolysis and oxidation of the thermo-char takes place. This reaction is exothermal with 
heat release H r The products are a-char (carbon) with density p ac , gas and black smoke with total density 

p^. Finally, combustion and oblation of the porous and friable a-char (carbon) layer formed at the 
previous step takes place at the TPS surface. This reaction is exothermal with heat release H [37]. The 

decomposition products influence essentially the base heating. Using earlier results [33] [36] the complex 
process of PU decomposition, formation of the thermo-char and a-char layer [37] can be reduced to the 
analysis of the following equations 
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Here R t , A t and A t are rates, coefficients, and activation energies of the corresponding reactions. 
Coefficients OC Tc , (3 ac show the fractions of solid thermo-char and a-char that are generated as a result of 
the pyrolysis of PU and thermo-char; (1 - a Tc ), (1 - / 3 ac ) are the fractions of the gas (see Eqs. (12)-(14)). 

The surface a-char layer has a porous structure and can oxidize (burn) in the atmosphere under a large 
enough heat flow [37], [39]. The rate of the surface processes and flux of the released gas are controlled 
by density of oxygen p Qx near the surface and can be presented as [37]: 

R s ,ac,i = expC-Zs / T s )p 0x (t), j,=a ef R saci , i = 0 2 ,C,C0,C0 2 (9) 

where p 0x (t ) is the local partial pressure of oxygen at the surface. Here the surface oxidation was described by the 
following brutto-reaction [37] 
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As a result of the PU pyrolysis a solid pore structure. The gas generated inside the pore structure is released due to 
pressure gradient at the surface. The decomposition process is very slow and the estimations show that the pressure 
drop is very small. In other words the gas pressure in the pores is close to the atmospheric pressure. As p atm <latm 

and T > T atm the gas density p < 1 kg / m 3 . Taking into account the fact that density and thermal conductivity of 

the gas are much smaller than for the solids we can neglect the gas contribution to heat balance that can be written 
as: 
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Mass and heat fluxes are determined for the multicomponent outgassing and smog are given by the following 
expressions 
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(13) 

(14) 


77 + 2 77 + 2 

where r s is the coordinate of moving surface of the burned a-char layer. The velocity of surface motion 
and mass flow of burned a-char are given by 
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Initial and boundary conditions: 
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T(x, 0 ) = T 0 , p p (x, 0 ) = p pu , p cT (x, 0 ) = 0 , p ac (x, 0 ) = 0 , j gi ( 0 ) = 0 ( 16 ) 
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3 

Typical values of the parameters used in the calculation are the following. PU parameters: p =(23 -f 30)kg/m , 
specific heat C = (1000 -f J/kg/K , thermal conduction 77 = (0.05 -fO. 75j W/m/K [34], [38], [39]. The 

p p 

10.1 4 5 8 

parameters of PU pyrolysis: A = 10 /sec, A =1.77x10 K, n=0.25, H = - 7.5x10 J/kg [36] or A = 1.686x10 /sec, 

pr p p pr 

A =1.62x10 K, n=0; H = - (2.93-3.47)xl0 J/kg [38]. Parameters of thermal char: p T =(l J rl0)kg/m , C t =1760 

J/kg/K, K = 0.05W/m/K . The parameters of the thermal char pyrolysis and generation of a-char layer: A = 

8.9 4 5 

10 /sec, n = - 0.08, E=l. 25x10 K; H t = 10 J/kg (with considering retardant) [36]. Heat of combustion of the a- 
char layer H burn =9.6xl0 6 J/kg [38]. It follows from the calculations that the a-char density is p =a ^ P c} p =12.8 

kg/m (a^fl ch =57%, p = 2 3 kg/m). 

Note that the ablation and combustion of the a-char layer takes place in two different scenarios at low and 
high altitudes. At low altitudes the combustions dominates and ablation is show. As a result the a-char 
layer is growing at the surface blocking heat conduction deep into the TPS. Therefore, the rate of PU 
pyrolysis and outgassing is small at low altitudes. An example of the results of calculations at low altitude 
is shown in the Figure 8. It can be seen from the figure that the front of the thermo-char paralysis 
propagates into the TPS with the rate v Tc ~ 0.125 mm/sec, while the recession rate of the TPS surface 
remains low v* ~ 0.046 mm/sec. 


Vetacirv of burning 



Figure 8: Kinetics of polyurethane foam pyrolysis for low altitudes. Radiation and convection heat flux Q 

4 2 3 

= 4-10 W/m , oxygen density p = 0.2kg/m . 
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Figure 9: Mass fluxes of gas, smog, oil, and char due to TPS ablation at high altitude. Weak burning and 

4 

strong convection ablation of a-char layer: was obtained for radiation and convection heat flux Q = 4-10 

2 3 

W/m and oxygen density p < 0.02kg/m . 


On the other hand, the higher is the altitude the lower is the concentration of the oxygen and the higher is 
velocity of the air trapped in the base. As a result, the ablation rate of the a-char layer [40] becomes large, 
giving rise to the increased flow of char particles with radius 1 to 10 pm and increased recession rate of 
the TPS surface to the values of the order of 1 mm/sec. This situation is illustrated in the Figure 9 where 
the mass fluxes of a number of components of the base bleeding are shown for TPS ablation at high 
altitude. 

The results of the calculations of the mass and heat fluxes to the base due to the TPS out-gassing can be 
used to include multi-component chemically active base bleeding into the CFD model of the base flow 
introduced in the previous sections. The results of the calculations of the base temperature as a function of 
the base bleeding parameters are presented and discussed below. 


7 Base overheating due to the base bleeding 

The analysis of the base environment in the presence of multi-component chemically reactive base 
bleeding was performed in a number of steps. First, the effect of outgassing with a given mass and heat 
fluxes was analyzed. Next, the char flow was added to the model. Finally, chemical reactions were 
included into the calculations. 


7.1 TPS outgassing effect on the base environment 

It was shown in the previous section that the TPS ablation and outgassing are complex functions of the 
altitude and base flow. In turn, the base flow depends strongly on the chemical composition, mass and 
heat fluxes of the TPS bleeding. To establish some of these dependences let us consider first the effect of 
outgassing on the base environment. Note that the effect of outgassing dominates at low altitudes where 
the recession rate of the char layer is small while the thermo-char pyrolysis inside the TPS is significant 
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Figure 10: (left) Relative changes of the base temperature (diamonds) as a function of the mass flow 
rate of the outgassing. (right) Variation of the relative base temperature as a function of the T bb . 
Straight lines are eye guiding interpolations of the numerical results. 


(see Sec. 6). The total gas flow rate and gas temperature estimated at low altitude using chemical kinetic 
model of section 6 are of the order of J bb = 0.015 kg/sec and T hb = 1200K. These estimations agree well 
with the experimental results reported in [32]. We now calculate the relative variations of the base 
temperature as a function of each one of these parameters while all other parameters of the model are kept 
constant. In particular, the char flow is assumed constant and equal to a small value J char = 1.5xl0" 4 
kg/sec. Examples of the corresponding dependences are shown in the Figure 10. It can be seen from the 
Figure 10 (left) that base temperature is approximately linear function of J bb in a wide range of variation 


of base bleeding mass flow rate. The linear dependence of T - (T b -T b0 )IAT on J bb is in agreement with 
predictions of the simplified model (7) of the temperature of the dead air in the trapped vortex in the base. 
Indeed, taking Q bb - J bb -c p • T bb , Q ch ~ 0, Q r = 0 and assuming that 
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one obtains linear dependence of the AT on J bb in a given time interval At. Similarly, the variation of the 
relative base temperature on the temperature of outgassing T bb is expected to be close to linear. This linearity is 

confirmed by the calculations of the T as a function of T bb shown in the Figure 10 (right). Note, however, 

that the relative error of these calculations is much larger than in the calculations of the dependence of T 
on J bb , because relative changes of the T bb are much smaller as compared to the changes of J bb . Note also 
that we have assumed above that T b0 is independent of the Q bb . This assumption is valid only for relatively 
small values of the Q bb . In general, the increase of the base temperature is expected to be a saturating 
function of J bb and T bb as will be discussed in details elsewhere. 


7.2 Char flow 


The effect of char flow on the base environment is even more complicated. On one hand, the heat 
exchange between the char particles and the environment is determined by a number of processes 
including the convective heat transfer, chemical reactions, and the absorption/emission of radiation at the 
particle surface. On the other hand, particles effectively block and absorb plume radiation, and may 
effectively escape from the base. In addition, the non-trivial characteristics of the particles behavior 
depend strongly on their radius. 

Consider, for example the rate of reaction at the particle surface given [28] 

R,=\n?,R, 
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Figure 11: (left) Variation of the relative base temperature as a function of the particle diameter D. Solid 
line is an eye guiding interpolation of the numerical results by function oc 1 ID. (right) Traces of the char 
particles colored by the temperature showing particles that escape the base flow. 


Where . is the rate of the particle surface depletion; A p is particle surface area; Yj is mass fraction of 


surface species rj is effectiveness factor; / R j is the rate of particle surface species reaction per unit area 
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where N = apparent order of reaction: 

K a ,=AT>e^‘ r ’ ) . 

It is clear that the rate of the surface reaction for the alpha-char is directly proportional to the total surface 
of the alpha particles. For constant mass flux ] C h = N4n/3-R 3 -p and reaction surface of the char particle A p 
- 4nR 2 the total reaction surface S = N-A p = 3J c iJRp and therefore the heat release are inversely 
proportional to the char radius. 

Similarly, radiation and the convective heat transfer are proportional to the particle surface area A p as 
shown in the equation (6). 

To verify the dependence of the convective heat exchange on the particle radius the variation of the 
relative base temperature was calculated for constant total char flow rate with different diameter of the 
char particles, while all other flow parameters were kept constant. The results of these calculations are 
shown in the Figure 11 (left). It can be seen from the figure that the measured decrease in the base 
temperature with an increase of the char diameter follows very closely the 1/R type of behavior. We have 
also verified the dependence of the particles escape rate from the base flow as a function of particle 
diameter. The results of this analysis demonstrate that the escape rate is larger for the particles of larger 
diameter. An example of the particles traces is shown in the Figure 11 (right). 

It is known from the experiment that the alpha-char particles radius can vary between 0.5 and 50 p. 
Therefore, it will be important in the future research to estimate distribution of alpha-char radius 
experimentally to estimate the effect of alpha-char burning on the base temperature. 

The results of the analysis discussed in this section show that the calculated effects of the base bleeding 
on the base temperature are consistent with theoretical expectations and with the experimental 
observations of the elevated base temperature during the test flight. In the next section we briefly 
summarize the effect of the chemical reactions and char burning on the base overheating. 


7.3 Chemical reactions and char burning in the base flow 

Multiple gas phase reactions occur in the base with the outcome depending on the operating conditions. 
Several variables seem to affect the gasification process, product composition, and distribution, including 
temperature, pressure, velocity, and oxygen mass fraction. These parameters are interrelated. The 
gasification rate is also temperature-dependent. The equilibrium points of the gasification reactions can be 
shifted by changing the temperature. The first-order rate constant of gasification was found to increase 
with temperature in accordance with the Arrhenius equation. It was assumed that the mixing of the 
components in the turbulent base flow was limiting the reaction rates. Accordingly, the eddy dissipation 
model was used to model chemical reaction in the gas phase with a number of limiters described in [28]. 
The simplified char combustion process includes three steps of devolatilization, char oxidation and gas- 
mixture oxidation among which the first two are heterogeneous (multi-phase) and the others are the 
homogeneous (single-phase) reaction. The eddy dissipation combustion model for homogeneous reactions 
was used in combination with a finite rate chemistry model for heterogeneous reactions at the particle’s 
surface. 



Figure 12: (left) Relative base temperature as a function of the model fidelity for two points on the flight 
trajectory: 1 - single-phase flow; 2 - A1 particles; 3 - radiation; 4 - char burning and chemical reactions 
(18% of particles with diameter D = 1 pm); 5 - char burning and chemical reactions (64% of 1 pm 
particles). Dashed lines indicate temperature measured during the flight, (right) Temperature profiles of 
the free stream plume interaction for multi-phase flow that includes A1 particles into multi-component 
plume flow with frozen chemistry and char particles with multi-component flow with homogeneous and 

heterogeneous reactions in the base flow. 

The results of the calculation of the base temperature that include chemical reactions in the gas phase and 
char combustion are shown in the Figure 12 (left) for two points on the flight trajectory. Red squares 
correspond to calculations performed for the point earlier in the flight with different levels of fidelity. In 
the figure notations 1 corresponds to the calculations performed for multi-component flow with a single 
gas phase after mesh adaptation. 2 corresponds to the calculations with A1 particles included into the 
plume flow. The last two members 4 and 5 correspond to the multi-component, multi-phase flow with 
chemical reactions included in the base flow both in gas phase and at the char particles surface (the plume 
flow assumes frozen chemistry). The difference between 4 and 5 is relative contributions of the char 
particles with diameter 1 pm. In case 4 particles with diameter D ~ 1 pm represent 18% of the total mass 
flow rate of the char flow. In case 5 they contribute 64% into the total mass flow rate. The dashed lines 
indicate the base temperature measured for given points of the flight trajectory during the test. It can be 


seen from the figure that the temperature calculated for the case of the base bleeding with for 
homogeneous and heterogeneous chemical reactions is approaching to the base temperature measured 
during the flight from below. 

The remaining gap between measured and calculated temperatures may have multiple origins. One of the 
contributions to the base temperature neglected so far was the plume radiation. The estimations show that 
the contribution of the radiation absorption into the elevated base temperature is of the order of 10K. 
Preliminary calculations of the model of the single-phase base flow that include radiation without 
chemical reactions are shown in Figure 12 (left) case 3. The results of these calculations are in agreement 
with the theoretical estimations. Another feature that may contribute to the elevated base temperature that 
was not included into the current model is an asymmetry of the base flow. Preliminary results of the 
calculations of the base flow in 90° sector that take into account flow asymmetry are shown in Figure 12 
(right) and will be reported elsewhere. The instabilities of the plume base flow interaction during e.g. 
maneuver during the flight may also result insignificant elevation of the base temperature. Finally, we 
note that the mass fluxes related to the TPS ablation are known only approximately and further research is 
required for more accurate predictions. Despite of remaining small discrepancies it can be concluded that 
the model of the base flow with chemically active base bleeding demonstrates correct trends in the base 
flow temperature that are consistent with observations and provide an explanation for one of the main 
sources of the heat flux into the base. 

8 Conclusion and Future Work 

Theoretical analysis and CFD simulations performed in this research revealed that the base heating 
anomalies observed in a recent test flight can be attributed to the following factors: small volume of the 
trapped base air in a single motor first stage configuration; significant base bleeding of hot chemically 
active components and char trapped in the base; and strong radial asymmetry of the base flow with small 
region of plume recirculation. An additional important factor that influence the overheating measured 
during the flight is a strongly non-uniform distribution of the hot flow in the base, which is mainly 
concentrated along the base walls for both base bleeding and plume recirculation. 

It was shown the base flow model with chemically active base bleeding can provide an explanation for 
one of the main sources of the heat flux into the base. The remaining sources of the heat fluxes that may 
contribute to significant elevation of the base temperature are plume radiation and transient instabilities of 
the plume flow. The development of the full 3D model of unsteady base flow that takes into account 
radiation is under way. The methodology of the base flow analysis developed in this research is currently 
being transferred, extended, and adopted for base flow in nozzle clusters. 

Nomenclature 

T = stress tensor 
ji - fluid viscosity 
77 = effectiveness factor 
K- thermal conductivity 
a =Stephan-Boltzman constant 
fi t - is the turbulent viscosity 
ju e ff =effective fluid viscosity 

R. = the rate of particle surface species reaction per unit area 

R- = the rate of the particle surface depletion 

A =rates coefficient 
A p = particle surface area 



B = non-dimensional pyrolysis gas rate at the surface 
B t =pre-exponential factor for the ith component 
C D = drag coefficient 
ch- Char 

Cp = particle concentration 
= solid material specific heat 
c v = specific heat at constant volume 
e = internal energy 
H= pyrolysis enthalpy 
h- specific enthalpy 
H c total enthalpy of the virgin material 
h re f - reference enthalpy 
H v total enthalpy of the charred material 
H w wall enthalpy 
J = mass flux of species 
k = turbulent kinetic energy 
m = mass of the particles 
Ma = Mach number 
N = apparent order of reaction 
rii = outward unit normal vector at surface i 
Nu =Nusselt number 
p = the static pressure. 
p' - total pressure 
Pr = Prandtl number 
q = radiative heat flux, W/m 1 2 3 4 
Q - the heat flux 
R - gas constant 
r- Reaction rate 
Re = Reynolds number 
R ex = radius of nozzle exit 
T - temperature 
t - time 

Ui - the ith Cartesian component of the velocity 
A = activation energies 
s = turbulent dissipation 
X = Conductivity 
p 7 = fluid density 
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